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Abstract 



The problem of Brownian motion in a periodic potential, under the 
influence of external forcing, which is either random or periodic in time, 
is studied in this paper. Multiscale techniques are used to derive general 
formulae for the steady state particle current and the effective diffusion 
tensor. These formulae are then applied to calculate the effective diffusion 
coefficient for a Brownian particle in a periodic potential driven simulta- 
neously by additive Gaussian white and colored noise. Our theoretical 
findings are supported by numerical simulations. 

1 Introduction 

Particle transport in spatially periodic, noi sy s ystems ha s att racted considerable 
attention over the last decades, see e.g. |33| Ch. 11], |3U] and the references 
therein. There are various physical systems where Brownian motion in periodic 
potentials plays a prominent role, such as Josephson junctions pQ, surface dif- 
fusion ITSlfSl and superionic conductors While the system of a Brownian 
particle in a periodic potential is kept away from equilibrium by an external, 
deterministic or random, force, detailed balance does not hold. Consequently, 
and in the absence of any spatial symmetry, a net particle current will appear, 
without any violation of the second law of thermodynamics. It was this funda- 
mental observation |21| that led to a revival of interest in the problem of particle 
transport in periodic potentials with broken spatial symmetry. These types of 
non- equilibrium systems, which are often called Brownian motors or ratchets, 
have found new and exciting applications e.g as the basis of theoretical mod- 
els for various intracellular transport processes such as molecular motors [S]- 
Furthermore, various experimental methods for particle separation have been 
suggested which are based on the theory of Brownian motors 0]. 

The long time behavior of a Brownian particle in a periodic potential is 
determined uniquely by the effective drift and the effective diffusion tensor 
which are defined, respectively, as 
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Here x(t) denotes the particle position, (•) denotes ensemble average and ® 
stands for the tensor product. Indeed, an argument based on the central limit 
theorem Ch. 3], ^B] implies that at long times the particle performs an 
effective Brownian motion which is a Gaussian process, and hence the first two 
moments are sufficient to determine the process uniquely. The main goal of 
all theoretical investigations of noisy, non-equilibrium particle transport is the 
calculation of and J5J. One wishes, in particular, to analyze the dependence 
of these two quantities on the various parameters of the problem, such as the 
friction coefficient, the temperature and the particle mass. 

Enormous theoretical effort has been put into the study of Brownian ratchets 
and, more generally, of Brownian particles in spatially periodic potentials 30 . 
The vast majority of all these theoretical investigations is concerned with the 
calculation of the effective drift for one dimensional models. This is not surpris- 
ing, since the theoretical tools that are currently available are not sufficient for 
the analytical treatment of the multi-dimensional problem. This is only possi- 
ble when the potential and/ or noise are such that the problem can be reduced 
to a one dimensional one |9] . For more general multi-dimensional problems one 
has to resort to numerical simulations. There are various applications, however, 
where the one dimensional analysis is inadequate. As an example we mention 
the technique for separation of macromolecules in microfabricated sieves that 
was proposed in 6 . In the two dimensional setting considered in this paper, 
an appropriately chosen driving force in the y direction produces a constant 
drift in the x direction, but with a zero net velocity in the y direction. On the 
other hand, a force in the x direction produces no drift in the y direction. The 
theoretical analysis of this problem requires new technical tools. 

Furthermore, the number of theoretical studies related to the calc ulation of 
the effective diffusion tensor has also been scarce [201 1311 1321 155] . In these 
papers, relatively simple potentials and/or forcing terms are considered, such 
as tilting periodic potentials or simple periodic in time forcing. It is widely 
recognized that the calculation of the effective diffusion coefficient is technically 
more demanding than that of the effective drift. Indeed, as we will show in this 
paper, it requires the solution of a Poisson equation, in addition to the solution 
of the stationary Fokker-Planck equation which is sufficient for the calculation of 
the effective drift. Diffusive, rather than directed, transport can be potentially 
extremely important in the design of experimental setups for particle selection 
|3()l Sec 5.11] [22]. It is therefore desirable to develop systematic tools for 
the calculation of the effective diffusion coefficient (or tensor, in the multi- 
dimensional setting) . 

From a mathematical point of view, non-equilibrium systems which are sub- 
ject to unbiased noise can be modelled as non-reversible Markov processes [2E] 
and can be expressed in terms of solutions to stochastic differential equations 
(SDEs). The SDEs which govern the motion of a Brownian particle in a periodic 
potential possess inherent length and time scales: those related to the spatial 
period of the potential and the temporal period (or correlation time) of the 
external driving force. From this point of view the calculation of the effective 
drift and the effective diffusion coefficient amounts to studying the behavior 
of solutions to the underlying SDEs at length and time scales which are much 
longer than the characteristic scales of the system. A systematic methodology 
for studying problems of this t ype, wh ich is based on scale separation, has been 
developed many years ago pfl 1241 l2"5"] . The techniques developed in the afore- 
mentioned references are appropriate for the asymptotic analysis of stochastic 
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systems (and Markov processes in particular) which are spatially and/or tempo- 
rally periodic. The purpose of this work is to apply these multiscale techniques 
to the study Brownian motors in arbitrary dimensions, with particular emphasis 
to the calculation of the effective diffusion tensor. 

The rest of this paper is organized as follows. In section [21 we introduce the 
model that we will study. In section 01 we obtain formulae for the effective drift 
and the effective diffusion tensor in the case where all external forces are Markov 
processes. In section0]we study the effective diffusion coefficient for a Brownian 
particle in a periodic potential driven simultaneously by additive Gaussian white 
and colored noise. Section[S]is reserved for conclusions. In Appendix A we derive 
formulae for the effective drift and the effective diffusion coefficient for the case 
where the Brownian particle is driven away from equilibrium by periodic in 
time external fluctuations. Finally, in appendix B we use the method developed 
in this paper to calculate the effective diffusion coefficient of an overdamped 
particle in a one dimensional tilted periodic potential. 



2 The Model 

We consider the overdamped d-dimensional stochastic dynamics for a state 
variable x(t) £ R d [301 sec. 3] 

yx(t) = -VV(x(t)J(t)) + y(t) + y/2jkjfZ(t), (3) 

where 7 is the friction coefficient, ks the Boltzmann constant and T denotes the 
temperature. £(t) stands for the standard d-dimensional white noise process, 

(&(*)> =0 and (&(*)&(«)> =*«<*(*-«). i,j=l,...d. 

We take f(t) and y{t) to be Markov processes with respective state spaces 
Ef, E y and generators £/, C y . The potential V(x, f) is periodic in x for every 
/, with period L in all spatial directions: 

V(x + LeiJ) = V(x,f), i=l,...,d, 

where {ei}f =1 denotes the standard basis of R d . We will use the notation Q = 

The processes f(t) and y(t) can be continuous in time diffusion processes 
which are co nst ructed as solutions of stochastic differential equations, dichoto- 
mous noise ^3 Ch. 9], more general Markov chains etc. The (easier) case 
where f(t) and y(t) are deterministic, periodic functions of time is treated in 
the appendix. 

For simplicity, we have assumed that the temperature in @ is constant. 
However, this assumption is with no loss of generality, since eqn. (|2J) with 
a time dependent temperature can be mapped to an equation with constant 
temperature and an appropriate effective potential |3(J1 sec. 6]. Thus, the above 
framework is general enough to encompass most of the models that have been 
studied in the literature, such as pulsating, tilting, or temperature ratchets. We 
remark that the state variable x(t) does not necessarily denote the position of 
a Brownian particle. We will, however, refer to x(t) as the particle position in 
the sequel. 

The process {x(t), f{t),y(t)} in the extended phase space M. d x Ef x E y is 
Markovian with generator 

C = F(x, /, y) ■ V x + DA X + C f +£ y , 
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where D := and 

7 

%/.l/) = -(-W(i,/) + y). 

7 

To this process we can associate the initial value problem for the backward 
Kolmogorov Equation |2.'il Ch. 8] 

— = Cu, u(x,y,f,t = 0)=u in (x,y,f). (4) 

which is, of course, the adjoint to the Fokker-Planck equation. Our derivation 
of formulae for the effective drift and the effective diffusion tensor is based on 
singular perturbation analysis of the initial value problem (TjJ . 

3 Multiscale Analysis 

In this section we derive formulae for the effective drift and the effective diffusion 
tensor for x(t), the solution of Let us outline the basic philosophy behind 
the derivation of formulae (|I3|I and (|I8|I . We are interested in the long time, 
large scale behavior of x(t). For the analysis that follows it is convenient to 
introduce a parameter e <C 1 which in effect is the ratio between the length 
scale defined through the period of the potential and a large "macroscopic" 
length scale at which the motion of the particle is governed by and effective 
Brownian motion. The limit e — > corresponds to the limit of infinite scale 
separation. The behavior of the system in this limit can be analysed using 
singular perturbation theory. 

We remark that the calculation of the effective drift and the effective diffu- 
sion tensor are performed seperately, because a different re-scaling is needed in 
each case. This is due to the fact that advection and diffusion have different 
characteristic time scales. 

3.1 Calculation of the Effective Drift 

The backward Kolmogorov equation reads 

du{x ^ Lt) = (F(x, /, y) ■ V x + DA X + C f + C y ) u(x, y, /, t). (5) 

We re-scale a space and time in Q according to 

x —>■ ex, t —y et 
and divide through by e to obtain 

% = 7 ( F /' y) ■ v * + eDA * + c f+ £ y) ut - 

We solve (0 pertubatively by looking for a solution in the form of a two-scale 
expansion 

u e {x, /, y, t) = u Q (x, ~,f,y,t)+ em (x, -,f,y,tj+ e 2 u 2 (x,-,f,y,tj+.... 

(7) 
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All terms in the expansion J7J are periodic functions of z — x/e. From the chain 
rule we have 

V X ^V X + -V Z . (8) 

Notice that we do not take the terms in the expansion JSJ) to depend explicitly 
on t/e. This is because the coefficients of the backward Kolmogorov equation 
do not depend explicitly on the fast time t/e. In the case where the fluctuations 
are periodic, rather than Markovian, in time, we will need to assume that the 
terms in the multiscale expansion for u e (x,t) depend explicitly on t/e. The 
details are presented in the appendix. 

We substitute now Q) into (JSJ, use JSJ and treat x and z as independent 
variables. Upon equating the coefficients of equal powers in e we obtain the 
following sequence of equations 

C u = 0, (9) 

£ "i = -Ciu a + -— , (10) 



where 

C =F{z,f,y)-\7 z + DA z + C y + C f (11) 

and 

A =F(z,f,y)- W X + 2DW Z W X . 

The operator £q is the generator of a Markov process on Q x E y x Ef. In order 
to proceed we need to assume that this process is ergodic: there exists a unique 
stationary solution of the Fokker-Planck equation 

£* oP (z >y> f) = 0, (12) 

with 

p(z,y, f)dzdydf 

QxEyXEf 

and 

C* P = V 2 • (F(z, f, y) P ) + DA zP + £* u p + £}p. 

In the above C*j and C* are the Fokker-Planck operators of / and y, respectively. 
The stationary density p{z, y, f) satisfies periodic boundary conditions in z and 
appropriate boundary conditions in / and y. We emphasize that the ergodicity 
of the "fast" process is necessary for the very existence of an effective drift and 
an effective diffusion coefficient, and it has b een tacitly assumed in all theoretical 
investigations concerning Brownian motors . 

Under the assumption that l|12|) has a unique solution eqn. implies, by 
Fredholm alternative, that uq is independent of the fast scales: 

mo = u(x, t). 

Eqn. 1)10(1 now becomes 

C oUl = - F(z, y, f) ■ V x u{x, t). 
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In order for this equation to be well posed it is necessary that the right hand 
side averages to with respect to the invariant distribution p(z, f, y). This leads 
to the following backward Liouville equation 

du(x.t) _ . . 

— ^ = U efr W x u{x,t), 

with the effective drift given by 

Ueff = / F(z, y, f)p(z, y, f) dzdydf 

J QxE y xEf 
1 



(-VV(x,f)+y)p(z,y,f) dzdydf. (13) 

7 JQxEy XEj 

3.2 Calculation of the Effective Diffusion Coefficient 

We assume for the moment that the effective drift vanishes, U e ff — 0. We 
perform a diffusive re-scaling in (J5J 

x — > ex, t — > e 2 t 

and divide through by e 2 to obtain 

1 (F(^,f,yj -V x + eDA x +£/ + £„) u\ (14) 



~dt 



We go through the same analysis as in the previous subsection to obtain the 
following sequence of equations. 

C u = 0, (15) 
C ux = -£iu , (16) 
£oU2 = —CiUi - C2U0, (17) 



where Cq and C\ were defined in the previous subsection and 

d 
at 

Equation Ijl5|l implies that uo — u(x,t). Now (|16(l becomes 

CqUi = -F(z, y, /) • V x u(x, t). 

Since we have assumed that U e ff = 0, the right hand side of the above equation 
belongs to the null space of Cq and this equation is well posed. Its solution is 

ui(x, z, f, y, t) = x(z, y, /) • V^u(x, t), 

where the auxiliary field x( z , y, f) satisfies the Poisson equation 

-C x(z,y,f) = F(z,y,f) 



G 



with periodic boundary conditions in z and appropriate boundary conditions in 
V and /. 

We proceed now with the analysis of equation i|17[l • The solvability condition 
for this equation reads 

(— C\Ui — C2U0) dzdydf = 0, 

XEyXEf 

from which, after some straightforward algebra, we obtain the limiting backward 
Kolmogorov equation for u(x, t) 

du{x,t) _ x - ^ eff d 2 u(x,t) 



J2 D 



dt ^-r 1 v dxidxj 

i,3 = l ■> 



The effective diffusion tensor is 



Df = D5 ij + (F i (z,yJ) X j (z,yJ)) p + 2D(p^^l^ , (18) 

where the notation (-) p for the averaging with respect to the invariant density 
has been introduced. 

The case where the effective drift does not vanish, Ueff 0, can be reduced 
to the situation analyzed in this subsection through a Galilean transformation 
with respect to U e f / 1 ■ The effective diffusion tensor is now given by 



D 



eff 



D5 l3 + ({Fi(z,y,f) - Ul ff )x 3 {z,yJ)) p 

+2P( ^% y,/) ) , (19) 

and the field x{ z i /> y) satisfies the Poisson equation 

-Lox = F(z >y> f)-U eff . (20) 

4 Effective Diffusion Coefficient for Correlation 
Ratchets 

In this section we consider the following model Ej 



7±(t) = -VV(x(t)) + y(t) + ^2 1 k B T^{t), (21a) 



m = --y(t)+J—at), (21b) 

T V T 



1 In other words, the process x^e)(t) := t (x(t/t 2 ) — e 2 U e fftj converges to a mean zero 
Gaussian process with effective diffusivity given by CHI 
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where £(t) and £(t) are mutually independent standard <i-dimensional white 
noise processes. The potential V(x) is assumed to be L-periodic in all spa- 
tial directions The process y(t) is the d-dimensional Onrstein-Uhlenbeck (OU) 
process which is a mean zero Gaussian process with correlation function 

_ |t-a[ 

{yi{t)Vj(s)) = Sijcre - , i,j = l,...,d. 

Let z(t) denote the restriction of x(t) to Q = [0, 2ir] d . The generator of the 
Markov process {z(t), y(t)} is 

C = -(-V z V(z) + y) • V z + DA Z + ^(-y-V y + aA y ) 

with D := £a£. Standard results from the ergodic theory of Markov processes 

see e.g. [3] ch. 3] ensure that the process {z(t),y(t)} 6 Q x M. d , with generator 
C is ergodic and that the unique invariant measure has a smooth density p(y, z) 
with re spect to the Lebesgue measure. This is true even at zero temperature 
Hence, the results of section El ap ply: the effective drift and effective 



diffusion tensor are given by formulae T ijT3l) and ltTH|) . respectively. Of course, 
in order to calculate these quantities we need to solve equations (|12J) and H20JI 
which take the form: 

--V 2 • {(-V z V(z) + y)p{y, z)) + DA zP (y, z) + - (V. y ■ {yp{y, zj) 
7 T 



+aA y p(y,z)) = 



and 



--(-V z V(z) + y) ■ W zX (y, z) - DA zX (y, z) 
7 

-\{-V V ?/X(y, z) + crAyx(y 7 z)) = I(_v,V(z) + y)-U. 

The effective diffusion tensor is positive definite. To prove this, let e be a unit 
vector in M. d , define f — F ■ e, u — U e ff ■ e and let <j> := e ■ x denote the unique 
solution of the scalar problem 

-C(j>=(F-U)-e=:f-u, 4>{y, z + L) = (f>(y, z), {<j>) p = 0. 

Let now h(y, z) be a sufficiently smooth function. Elementary computations 
yield 

C*{hp) = -pCh + 2DV Z ■ (pV„h) + • (pV y h) . 

We use the above calculation in the formula for the effective diffusion tensor, 
together with an integration by parts and the fact that (<p(y, z)) p = 0, to obtain 

e-D efr e = D + {f<j)) p + D{e-V z 4>) p 

= D+{u<j>) p -(cj)£cj ) ) p + 2D{e-V y cb) p 

= D + D{\V Z ^\ 2 ) P + 2D{e ■ V,» p + -(|V^| 2 ) P 

T 

= ^(| e + V z 0| 2 ) p + -(|V^| 2 ) p . 

T 



From the above formula we see that the effective diffusion tensor is non-negative 
definite and that it is well defined even at zero temperature: 



e-D eff (T = 0).e=-(\V v <f>(T = 0)\ 2 ) p . 

Although we cannot solve these equations in closed form, it is possible to cal- 
culate the small r expansion of the effective drift and the effective diffusion 
coefficient, at least in one di mension . Indeed, a tedious calculation using singu- 
lar perturbation theory, e.g. ^3 |2J] yields 

Ueff = 0(r% (22) 

and 

Drf/ -^( D+ „( 1 + ^(|-|))) +0(A <*, 

In writing eqn. (|23[1 we have used the following notation 

f L _vm ~ f L vw 

Z = I e o dz, Z = / e D dz, 
J o Jo 

Zi= [ V(z)e~ Z ^ 1 dz, Z 2 = [ V(z)e X ^ 1 dz. 
Jo Jo 

It is relatively straightforward to obtain the next order correction to i|23[l : the 
resulting formula is, however, too complicated to be of much use. 

The small r asymptotics for the effective drift were also studied in [21 El for 
the model considered in this section and in 7,36 when the external fluctuations 
are given by a continuous time Markov chain. It was shown in [71 13fi| that, for 
the case of dichotomous noise, the small r expansion for U e /f is valid only for 
sufficiently smooth potentials. Indeed, the first non-zero term-of order C(t 3 )- 
involves the second derivative of the potential. Non-smooth potentials lead to 

an effective drift which is O(r^). On the contrary, eqn. I|23[l does not involve 
any derivatives of the potential and, hence, is well defined even for non-smooth 
potentials. On the other hand, the C(r 2 ) term involves third order derivatives 
of the potential and can be defined only when V(x) £ C 3 (0, L). 

We also remark that the expansion (12311 is only valid for positive tempera- 
tures. The problem becomes substantially more complicated at zero tempera- 
ture because the generator of the Markov process becomes a degenerate differ- 
ential operator at T = 0. 

Naturally, in the limit as r — > the effective diffusion coefficient converges 
to its value for y = : 

ft„ = ff. (24) 

This is the effective diffusion coefficient for a Brownian p article m oving in a 
periodic potential, in the absence of external fluctuations |19l 13*7] . It is well 
known, and easy to prove, that the effective diffusion coefficient given by l|24|l 
is bounded from above by D. This not the case for the effective diffusivity of 
the correlation ratchet 12 1|) . 

We compare now the small r asymptotics for the effective diffusion coeffi- 
cient with Monte Carlo simulations. The results presented in figures 1 and 2 
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Figure 1: Effective diffusivity for 1|21|) with V(x) = cos(x) as a function of r, for 
a = 1, D = %p = 1, 7 = 1. Solid line: Results from Monte Carlo simulations. 
Dashed line: Results from formula (|23|l . 

were obtained from the numerical solution of equations l|21(l using the Euler- 
Marayama method, for the cosine potential V(x) = cos(ir). The integration 
step that was used was At — 10~ 4 and the total number of integration steps 
was 10 7 . The effective diffusion coefficient was calculated by ensemble averag- 
ing over 2000 particle trajectories which were initially uniformly distributed on 
[0,2tt]. 

In figure 1 we present the effective diffusion coefficient as a function of the 
correlation time r of the OU process. We also plot the results of the small 
r asymptotics. The agreement between theoretical predictions and numerical 
results is quite satisfactory, for r <C 1. We also observe that the effective 
diffusivity is an increasing function of r. 

In figure 2 we plot the effective diffusivity as a function of the noise strength 
a of the OU process. As expected, the effective diffusivity is an increasing 
function of a. The agreement between the theoretical predictions from (|23|l and 
the numerical experiments is excellent. 



5 Summary and Conclusions 

The problem of Brownian motion in a periodic potential and driven by external 
forces was studied in this paper. Using multiscale techniques we derived for- 
mulae for the effective drift and the effective diffusion tensor. We used these 
formulae to calculate the effective diffusion coefficient for a Brownian particle in 
a one dimensional periodic potential which is driven by additive colored noise. 
The theoretical predictions are in very good agreement with the results obtained 
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Figure 2: Effective diffusivity for l|21|) with V(x) = cos(x) as a function of 
a, for t = 0.1, D = = 1, 7 = 1. Solid line: Results from Monte Carlo 
simulations. Dashed line: Results from formula i|23|) . 



from direct numerical calculations. 

In order to streamline the presentation, we studied the cases of Markovian 
and time-periodic external fluctuations separately. It is straightforward, of 
course, to combine the two and obtain formulae for the effective drift and the 
effective diffusion coefficient when some of the processes are random and some 
periodic. Furthermore, the case where the forcing terms depend on the particle 
position can be treated easily within the framework developed in this paper. In 
order for our method to be applicable in this case it is necessary for the fast 
process to be ergodic for all values of the state variable x. We believe, thus, that 
our multiscale approach is general enough to encompass essentially all models for 
periodic potentials and external fluctuations that have been considered within 
the framework of Brownian motors. 

The multiscale technique employed in this paper provides us with a system- 
atic method for analyzing the long time behavior of solutions to the stochas- 
tic equations of motion. The assumption of scale separation upon which the 
method is based is always satisfied for systems with a finite number of charac- 
teristic length and time scales. Naturally, the results of the multiscale analysis 
are valid at length and time scales large compared to the largest length and time 
scal es of the system. For example, the time required for the correlation ratchet 
<|21l) to reach its asymptotic Brownian regime increases as the correlation time 
r increases. However, for finite r the system will always eventually perform an 
effective Brownian motion. 

In this paper we only presented formal calculations. We emphasize however 
that all the results presented in this paper can be justified rigorously, and in 
a very general setting. The key a ssu mptions are the ergodicity of the fast 
process generated by £q, defined in (|11|) . and the well-posedness of the Poisson 
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equation (|2UI) . Ergodicity implies the existence of the effective drift: this is a 
consequence of the stro ng law of large numbers. Ergodicity, together with well- 
posedness of eqn. (|2U[1 imply the existence of the effective diffusion coefficient: 
this is a consequence of the martingale central limit theorem, see e.g. ch. 
7]. Needless to say, checking the ergodicity of the fast process and the well 
posedness o f H2()|l are not trivial tasks. Results of this form are proved, for 
example, in |14|. 

In the case of stochastic external fluctuations, the calculation of the effective 
diffusion tensor requires the numerical solution of equations l|12l) and (|20|l . It 
is conceivable that this task might be as computationally demanding as direct 
Monte Carlo simulations. From this point of view the main advantage of our 
approach is that it provides a natural starting point for careful asymptotic 
analysis and rigorous study of the dependence of the effective drift and the 
effective diffusion coefficient on the various parameters of the problem. 

On the contrary, when the coefficients of eqn. are per iodi c in t ime , 
the calculation of D e ff requires the numerical solution of eqns. (|30[1 and (|37|) . 
These are PDEs with periodic coefficients and can be routinely solved with e.g. 
a spectral method. In this setting, our approach offers a clear computational 
advantage over Monte Carlo simulations. We also mention in passing that the 
various formulae for the effective diffusion coefficient that have be en derived 
in the literature ^1 El EH1 OHI can be obtained from equation l!3fi|) : they 
correspond to cases where equations (l3*TJ) and (|3"?|) can be solved analytically. An 
example-the calculation of the effective diffusion coefficient of an overdamped 
Brownian particle in a tilted periodic potential-is presented in appendix. Similar 
calculations yield analytical expressions for all other exactly solvable models that 
have been considered in the literature. 

In this paper we have studied the overdamped limit. There are various 
problems, however, where inertial effects are important and cannot be neglected. 
As exam ples we mention rocking ratchets in SQUIDS, [301 sec. 5.10], surface 
diffusion 18, 34], inertial particles |2 71 1 1 4j . Our framework can be generalized to 
include inertial effects. It is very interesting, then, to understand the dependence 
of the effective diffusion coefficient on the particle inertia. This problem is 
currently under investigation. 

Acknowledgements The author is indebted to Professors OR. Doering 
and A.M. Stuart for useful discussions and suggestions. 



A Time— Periodic Coefficients 

In this appendix we derive formulae for the mean drift and the effective diffusion 
coefficient for a Brownian particle which moves according to 

7 ±(t) = -W(x(t), t) + y(t) + y/2-yk B T{x{t),t)Z(t), (25) 

for space-time periodic potential V(x, t) and temperature T(x,t) > 0, and 
periodic in time force y(t). We take the spatial period to be L in all directions 
and the temporal period of V(x, t), T(x, t) and y(t) to be T. We use the notation 

Q = [0, L] d . Equation (1251) is interpreted in the Ito sense. This is the correct 
interpretation when ean. l|25)l is obt ained fr om the full phase space dynamics, 
in the limit as the inertia tend s to |17l I2(ij . Since the analysis is very similar 
to the one presented in |37] we will be brief. 
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A.l Calculation of the Effective Drift 

The backward Kolmogorov equation corresponding to 1(25(1 reads 

= F(x, t) ■ V x u{x, t) + kBTM A x u(x, t), (26) 
at 7 

where F(x,t) = ^(— W(x, y) + y(t)). We re-scale space and time in iffifl) 

according to x =— ► err, t — > et. We look for a solution of the resulting equation 
of the form 

u e (x, t) = u (x, —,t, + eui ^x, -,t,-^j + e 2 u 2 (x, — ,t, + . . . . (27) 

We treat x, z — x/e and t, r = t/e as independent variables. We emphasize 
that Uj(x, z, t, t) , i = 0, 1 , 2, . . . are periodic functions in both z and r. Upon 
substituting 127(1 into 1(26(1 , using the chain rule and equating the coefficients of 
equal powers in e we obtain the following sequence of equations 

Tl o u = °> (28) 
TZo^ = -Tlmo, (29) 



where 



and 



Uo = ^--C , Co = F(z,r)-V z + ^^A z 
or 7 

d , . 2k B T(z,T) 

n x = — -L u Ci = F (z,t -v a + - v ' v z v x . 
at 7 



The equation 

dp(z,r) i ^ _ w _ _^ , k B T(z,r 

_ 

with 



ft*p(z, r) = ^LLL + V z (F(z, r)p(z, t)) + ) ' J A z p(z, r) = 0, (30) 



r 

p(z, t) dzdr = 1 

and periodic boundary conditions in both z and r has a unique solution under 
appropriate regularity ass umptions on F{z,t), T(z,t) Thm. 3.10.1]. By 
Fredholm alternative, ((28(1 implies that 



uq = u(x, t). 

Eqn. 1(29(1 now becomes 

Tl oUl = _ T ) . v , u (x, t). 

at 

The solvability condition for this equation leads to the following backward Li- 
ouville equation 

du(x,t) . . 

— ^ — = ^e// ' v x u(x,t), 

with 



U eff = [ [ F(z,T)p(z,T)dzdT = 0. (31) 
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A. 2 Calculation of the Effective Diffusion Coefficient 

Without loss of generality we assume that the effective drift given by 
vanishes, U e ff = 0. We perform the diffusive re-scaling in equation 
x — > ex, t — » e 2 t. We go through the same analysis as in the previous sub- 
section to obtain the following sequence of equations. 

n u = 0, (32) 
IZqui = —Ciuq, (33) 
TZqU2 = -C%ui - C 2 u (34) 

where IZo, Co, Ci were defined in the previous subsection and 

d k B T(z,r) A 
C2 = dl 7 Ax - 

Equation l|32(l gives that uo = u(x,t). Equation (|33|l becomes 

U ui = —F(z, t) ■ V x u(x, t). 

Since we have assumed that U e ff = 0, this equation is well posed. Its solution 
is 

ui(x,z,t,r) = x(z,t) -V x u(x,t), 

where 

-TZ oX (z,t) = F(z,t). (35) 
The solvability condition for equation Ij34(l reads 



/ / (— C\U\ — C2U0) dzdr = 0, 
Jo Jq 

from which, after some algebra we obtain an evolution equation for u{x, t) 



ii 



du(x,t) d u(x,t) 

dt y dxidxj 



The effective diffusion coefficient is 

Aj = hi — (T(»,T)>,+ (J',(j,T)xj(z,T)) ) , 

7 \ d Vj I p 

with 

T r 

p(z, t) dzdr. 

When U e ff ^ 0, eqns. and become, respectively, 

-Kox{z,t) = F{z,r) - U e/f (37) 
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and 

Dii = S^{T{z,r)) p + {{Fi{z,r)-Ui Jf )^{z,r)) p 

B Effective Diffusion Coefficient for Tilted Pe- 
riodic Potentials 

In this appendix we use our method to obtain a formula for the effective diffusion 
coefficient of an overdamped particle moving in a one dimen sional t ilted periodic 
potential. This formula was first derived and analyzed in |32l I31| without any 
appeal to multiscale analysis. The equation of motion is 

x = -V'(x)+F + V2D£, (39) 

where V(x) is a smooth periodic function with period L, F and D > constants 
and £(t) standard white noise in one dimension. To simplify the notation we 
have set 7 = 1. 

The stationary Fokker-Planck equation corresponding to (|39ll is 

d x ((V'{x) ~ F) p(x) + Dd xP {x)) = 0, (40) 

with periodic boundary conditions. Formula l|13|) for the effective drift now 
becomes ^ 

U eff = I (-V'(x)+F)p(x)dx. (41) 
Jo 

The solution of eqn. gDJ is 29, Ch. 9] 

P{x) = \ dyZ + (y)Z_(x), (42) 

with 
and 

pL px-\-L 

Z= dx dyZ + (y)Z^(x). (43) 

JO Jx 

Upon using g5J> in iJUJl we obtain g§l Ch. 9] 

U eff = ^(l-e-^). (44) 

Our goal now is to calculate the effective diffusion coefficient. For this we first 
need to solve the Poisson equation l|2UI) which now becomes 

C X (x) := Dd xxX (x) + (-V'(x) + F)d xX = V'(x) -F + U eff , (45) 
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with periodic boundary conditions. Then we need to evaluate the integrals in 



pL pL 

D eff = D+ (~V'{x)+F-U eff )p(x)dx + 2D d xX {x)p{x) 
Jo Jo 



dx. 



It will be more convenient for the subsequent calculation to rewrite the above 
formula for the effective diffusion coefficient in a different form. The fact that 
p(x) solves the stationary Fokker-Planck equation, together with elementary 
integrations by parts yield that, for all sufficiently smooth periodic functions 

<f>(x), 

pL pL 

/ (j>(x)(—C(f>(x))p(x)dx=D / {d x <t>{x)fp{x)dx. 
Jo Jo 

Now we have 

D eff = D+ {-V'(x)+F-U eff )x{x)p(x)dx + 2D d xX {x)p{x)dx 
Jo Jo 

pL pL 

= D+ (~C X (x)) X (x)p(x)dx + 2D d xX (x)p(x)dx 
Jo Jo 

pL pL 

= D + D {d xX {x)f p{x) dx + 2D d xX (x)p(x) dx 
Jo Jo 

r L 

= D / (l + d xX (x)f p{x)dx. (46) 
Jo 

Now we solve the Poisson equation 11451 with periodic boundary conditions. We 
multiply the equation by Z-(x) and divide through by D to rewrite it in the 
form 

d x (d xX (x)Z^(x)) = -d x Z_{x) + ^Z_(x). 

We integrate this equation fr om x — L to x and use the periodicity of X (x) and 
V(x) together with formula (|44|) to obtain 



d xX (x)Z_(x)(l-e-^ ^-Z_(a;)(l-e-^)+|(l-e-^) J ^Z_(y)dy, 
from which we immediately get 



8 xX {x) + 1 = \ I Z_(y)Z + (x)dy. 



Z 



—L 



Sub stituting this into (|46|l and using the formula for the invariant distribution 
1)421) we finally obtain 

Deff = §>f( T +(x)) 2 I-(x)dx, (47) 

with 

rx px+L 

I +( x ) = Z-(y)Z + (x)dy and = / Z + {y)Z-(x) dy. 

J x — L J x 

Formula (|47|l for the effective diffusion coefficient (formula (22) in 31 ) is the 
main result of this appendix. 
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